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ABSTRACT 

Rapid UBVRI photometry of the ultracompact low-mass X-ray binary (LMXB) 
pulsar 4U 1626-67/KZ TrA has detected 130.4 mHz (7.67 s) optical pulsations 
in all hve bands. The optical pulsations, which are at the same freqnency as the 
X-ray pnlsations cansed by rotation of the highly-magnetized accreting nentron 
star primary, are nnderstood as reprocessing of the pnlsed X-ray emission in the 
accretion disk or on the snrface of the secondary. The optical pnlsed fraction is 
ronghly 6%, independent of wavelength, indicating that the optical emission is 
dominated by X-ray reprocessing. 

A weaker (1.5%) sideband, downshifted 0.395(15) mHz from the main optical 
pnlsation, is also present. This is consistent with a previonsly reported sideband 
which was downshifted 0.4011(21) mHz from the main pnlsation, corroborating 
the 42-min binary period proposed by Middleditch et al. (1981, ApJ, 244, 1001). 

A 0.048 Hz optical qnasi-periodic oscillation (QPO), corresponding to a previ¬ 
ously reported X-ray feature, was also detected in some of the observations, with 
a fractional RMS amplitnde of 3-5%. This is the first measnrement of an optical 
QPO in an X-ray binary pnlsar. 

I discnss constraints on the natnre of the mass donor and show that mass 
transfer via a radiatively-driven wind is inconsistent with the optical data. I also 
review the basic theory of X-ray-heated accretion disks and show that snch models 
provide a good fit to the optical photometry. If the effective X-ray albedo of 
LMXB accretion disks is as high as recently reported (p^ 0.9), then the optical 

data imply a distance of ~8 kpc and an X-ray Inminosity of ~ 10^^ erg s“^. 
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1. INTRODUCTION 

The optical emission from low mass X-ray binaries (LMXBs) is generally dominated by 
the reprocessing of X-rays in the accretion disk and/or the mass donor (see van Paradijs & 
McClintock 1995 for a review). The time history of the optical emission is thus a convolution 
of the X-ray intensity history with a function representing the spatial distribution of matter 
in the system. The most spectacular examples of this phenomenon are the optical novae 
accompanying soft X-ray transients, which can cause the optical counterpart to brighten by 
many magnitudes. 

Of greater interest are those systems whose X-ray intensity histories are more regularly 
or sharply modulated. Many LMXBs show periodic optical variability due to orbitally- 
modulated viewing of the X-ray heated mass donor, providing constraints on the binary 
inclination and the distribution of reprocessing material in the system (see van Paradijs 
& McClintock 1995 and references therein). Quasi-periodic oscillations (QPOs) at various 
frequencies have been observed in both the X-ray and optical emission from the candidate 
black hole binary GX 339-4 (Motch, Ilovaisky, & Chevalier 1982; Motch et ah 1983; Imamura 
et al. 1990). Correlated X-ray and optical bursts have been observed from several systems, 
and the time delay and smearing of the optical burst prohles with respect to the X-ray 
bursts have been used to probe the location and distribution of matter near the neutron star 
(Lewin, van Paradijs, & Taam 1993). Coherent optical pulsations have been detected from 
the LMXB pulsars Her X-1 (Davidsen et al. 1972; Middleditch & Nelson 1981), 4U 1626-67 
(Ilovaisky, Motch, & Chevalier 1978; Middleditch et al. 1981), and GX 1-1-4 (Jablonski et al. 
1997), and likewise provide an important probe of the binary parameters and the distribution 
of reprocessing material in the system. 

The LMXB 4U 1626-67/KZ TrA consists of a 7.66 s pulsar accreting from an extremely 
low-mass companion in an ultracompact binary with a separation of < 1 It-sec (Middleditch 
et al. 1981; Levine et al. 1988). The 7.66 s X-ray pulsations arise due to anisotropic 
accretion of matter on the surface of a rotating, highly-magnetized neutron star whose spin 
and magnetic dipole axes are misaligned. The spin frequency of the pulsar evolves on short 
(|z//z>| ^ 6000 yr) time scales due to accretion torques (Chakrabarty et al. 1997). The optical 
counterpart has a strong ultraviolet excess (McClintock et al. 1977), and optical pulsations 
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are detected with the same frequency as the X-ray pulsations (Ilovaisky et ah 1978). The 
optical emission is understood as primarily due to reprocessing of the incident X-ray flux 
by material in the accretion disk and/or the binary companion’s surface (Chester 1979; 
McClintock et ah 1980). The system shows strong, correlated X-ray/optical flares every 
~ 1000 s that are of undetermined origin (Joss, Avni, & Rappaport 1978; McClintock et al. 
1980; Li et al. 1980). A ~ 0.04 Hz QPO has been detected in the X-ray emission (Shinoda 
et al. 1990; Angelini et al. 1995), and probably arises through an interaction between the 
pulsar’s magnetosphere and the inner edge of the accretion disk (Alpar & Shaham 1985; 
Lamb et al. 1985). 

High-sensitivity timing of the optical pulsations detected an additional weak, persistent 
pulsation in a lower-frequency sideband of the “direct” (X-ray) pulse frequency (Middleditch 
et al. 1981). This sideband was attributed to X-ray reprocessing on the surface of the 
companion, so that the 0.4 mHz downshift is due to the different effective pulsar rotation 
rate observed in the companion frame due to the binary orbit|. The observed frequency 
shift implies a 42-min prograde orbit. X-ray timing measurements have repeatedly failed 
to detect periodic pulse arrival time delays or pulse frequency shifts due to a binary orbit, 
yielding an upper limit of a^isini ^ 10 It-ms for the projected radius of the neutron star 
orbit (Levine et al. 1988; Shinoda et al. 1990). For a 42-min orbit, the X-ray timing limits 
imply a mass function fx{M) < 10“®Mq, one of the smallest known for any stellar binary. 

This paper describes a program to conhrm the persistent orbital sideband reported by 
Middleditch et al. (1981) and also provides the hrst detailed study of the optical wavelength 
dependence of X-ray reprocessing in this system. As part of this program, the hrst example 
of an optical QPO in an X-ray pulsar system was discovered. A preliminary account of this 
work has appeared previously (Chakrabarty 1996). 


2. OBSERVATIONS 

High speed multicolor optical photometry of 4U 1626-67/KZ TrA was obtained on UT 
1995 May 26-28 (MJD 49863-49865) with the 4-m Blanco telescope at the Cerro Tololo Inter- 
American Observatory (CTIO) in La Serena, Chile. The observations used the Automated 
Single Channel Aperture Photometer (ASCAP) and a dry-ice-cooled Varian VPM-159A 
photomultiplier tube at the f/7.5 Ritchey-Chretien focus, along with the UB VRI fUter set 
described by Graham (1982). The Varian phototube, which has an InGaAsP photocathode, 
is sensitive over a wide wavelength range, with a quantum efficiency of 15% at 4000 A and 
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5% at 9000 A. All the observations were made through a 6.6 arcsec circular aperture. The 
data were recorded at 1 ms resolution and rebinned to 100 ms resolution for timing analysis. 
A log of the long timing observations of 4U 1626-67 is given in Table 1. Typical target and 
sky count rates are given in Table 2. The Graham (1982) photometric standard star E7-s 
was also observed at various airmasses on each night as a calibration source. 

Timing stability was maintained using a Kinemetrics Truetime 468-DC GOES satellite 
synchronized clock. This clock was synchronized using NIST timing signals transmitted via 
the NOAA Geostationary Operational Environmental Satellite 7 (GOES/West), located in 
geostationary orbit at 135° west longitude (Beehler & Lombardi 1990). The signals were 
corrected for the 267 ms mean path delay from the NOAA ground station at Wallops Island, 
Virginia up to GOES/West and down to the GTIO receiver (corresponding to a setting of 
57 ms for the propagation delay switch on the clock). According to the manufacturer spec- 
ihcations, the clock in this conhguration maintains absolute time accuracy within ±1.5 ms 
of UTG(NIST). However, an independent check of this accuracy (e.g., a cross-check with a 
Global Positioning System timing signal or optical pulse timing of a well-established astro- 
physical time standard like the Grab pulsar) was not available during these observations. 

Some of the long timing observations contained a strong instrumental signal at 60 Hz 
and its higher harmonics in the 1 ms data (see Table 1). Whenever the 60 Hz signal was 
present, it was accompanied by a second strong signal at 1.75 Hz and its higher harmonics, 
and the two signals beat against each other. These high frequency contaminants appeared 
intermittently in observations of several different targets. However, their presence did not 
affect any of the science analysis for 4U 1626-67, all of which was conhned to pulse frequencies 
below 1 Hz. 


3. RESULTS 

Absolute photometry of 4U 1626-67 on May 26 yielded the following magnitude and 
colors: V = 18.68(15), U - B = -1.20(17), B - V = 0.02(19), and V - R = 0.00(21). 
The individual magnitudes are given in Table 2. The /-band sky background observation 
was too short to allow determination of a reliable source magnitude, but the presence of 
7.66-s /-band pulsations, in phase with those in the bluer bands (see below), conhrms that 
the source was detected in the /-band as well. Assuming that the /-band pulsed fraction 
is also 6% (as in the other optical bands; see below), the source has / = 17.8(3). Within 
the uncertainties, interstellar reddening can be neglected for this high latitude source, since 
the total Galactic reddening in the direction of 4U 1626-67 is E{B — U) ~ 0.09 (Burstein & 
Heiles 1982). The UBV magnitudes reported here are virtually identical to those measured 



5 


on 1977 June 12 by McClintock et al. (1977), and they agree within 2a with those measured 
on 1977 August 17 by Grindlay (1978). However, the present i?-band measurement is 1.8(4) 
magnitudes brighter than previously reported by Grindlay (1978). 

The observed time series were transformed to the solar system barycenter frame using 
the Jet Propulsion Laboratory DE-200 solar system ephemeris (Standish et al. 1992). Power 
spectra of the individual observations in all hve bands showed strong 7.66 s pulsations, with 
a mean barycentric pulse frequency of 130.4398(7) mHz. This is consistent with the X-ray 
pulsation frequency measured at the midpoint of the GTIO observations by the BATSE 
all-sky monitor on the Compton Gamma Ray Observatory (Ghakrabarty et al. 1997). To 
facilitate a study of the pulse shape and phase relationship of the optical pulsations at dif¬ 
ferent wavelengths, a precise pulse phase ephemeris was derived from a subset of the BATSE 
observations (1995 May 9 - 1995 June 11) contemporaneous with the GTIO observations. 
Based on these data, the barycentric X-ray pulse phase during the GTIO observations is well 
£t by 

0(t) = 00 + z^o(t - to) + ^z>(t - ^o)^ (1) 

where 0o = 0.936758, uq = 0.1304376013765 Hz, z> = -6.32044 x lO'^^ Hz s-0 R = 1995 May 
27.5 TDB, and pulse phase is defined relative to the intensity minimum of the fundamental 
harmonic of the 20-60 keV BATSE pulse profile. 

The GTIO timing observations for each of the hve filters were folded according to the 
timing model in Equation (1). The resulting pulse profiles, along with the BATSE 20-60 
keV pulse profile, are shown in Figure 1 relative to the BATSE pulse phase. The mean 
optical pulse and the BATSE hard X-ray pulse agree in phase to within 5% (i.e., a relative 
time delay of ^ 0.4 s), consistent with previous measurements (McGlintock et al. 1980). 
Both the BATSE hard X-ray pulse and the optical pulse shapes have a dominant sinusoidal 
component, although the optical pulses clearly have additional harmonic content as well. 
The optical pulse shapes do not resemble any of the X-ray pulse shapes in the 2-13 keV 
range, all of which have an unusual and distinctive close double-peak structure (Levine et al. 
1988). Instead, the optical pulse shapes resemble the low-energy (1-2 keV) and high-energy 
(15-60 keV) pulse shapes. This is consistent with previous observations by McGlintock et al. 
(1980), who concluded that the 1-2 keV X-rays were chiehy responsible for the reprocessed 
optical emission. The pulsed fraction, independent of wavelength, is approximately 6% (see 
Table 2), suggesting that the optical emission is dominated by X-ray reprocessing. 

Figure 2 shows the Fourier power spectrum of a 2-hour 7/-band observation (run 15U). 
The spectrum has been rebinned into uniform logarithmic frequency intervals, the mean 
noise level due to Poisson counting statistics subtracted off, and the power level normalized 
relative to the double-sided mean source power (e.g., Miyamoto et al. 1994). In addition 
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to the 130 mHz pulsation, a 0.048 Hz QPO was also detected. This feature was £t with a 
Lorentzian prohle (e.g., Shibazaki & Lamb 1987), 


Pq¥o{v) CC 




TT 


( 2 ) 


(Az//2)2 + (z/ - z/QPo)^’ 
where z^qpo is the centroid frequency of the QPO and Az/ is its full width at half-maximum 
(FWHM). The fractional root-mean-squared (RMS) amplitude of the QPO was computed by 
integrating equation (|^) over the FWHM of the feature and taking the square root, yielding 
3.4% for this observation. The QPO is also detected in the B and R-band observations. 
The best-fit parameters for the QPO features are included in Table 2. These measurements 
provide the hrst detection of an optical QPO from an X-ray pulsar. A similar feature is 
observed at X-ray wavelengths (Shinoda et ah 1990; Angelini et ah 1995; Chakrabarty et 
ah 1997, in preparation). For comparison, the Fourier power spectrum of the 2-6 keV data 
from a 71 ks observation of 4U 1626-67 with the Proportional Counter Array on the Rossi 
X-Ray Timing Explorer (RXTE) is shown in Figure 3. This observation was made on 1996 
February 11. The QPO is clearly visible in the X-ray data with an RMS strength of 6%, 
along with several harmonics of the coherent 7.66 s pulsation. (The higher harmonics are 
very strong due to the very sharp 2-6 keV X-ray pulse shape of 4U 1626-67.) The RXTE 
data will be discussed in further detail elsewhere (Chakrabarty et ah 1997, in preparation). 


To study the detailed structure of the optical timing properties close to the pulsation 
frequency, an overresolved power spectrum was constructed by computing the Fourier power 
Pk+e = \Hk+e\^ at fractional bin intervals according to 

Af-l N-l 

Hk+e = E = E (3) 

i=o i=o 

with 0 < e < 1, where y^^'^ is just the original time series hj multiplied by complex phase 
factors . Note that while the N shifted power spectrum bins {Pk+e] are statistically 

independent of each other, they are not independent of the original power spectrum {Pk} 
(although the covariances are easily computed; see, e.g., Jenkins & Watts 1968). Generaliz¬ 
ing this process, one can construct an n-times overresolved power spectrum by interleaving 
the bins from the original power spectrum and rz, — 1 frequency-shifted power spectra. Op¬ 
erationally, this is equivalent to padding the original A-point zero-mean time series with 
{n — 1)N zeros and computing the power spectrum of the new rzA-point padded time series. 
Use of an oversampled power spectrum gives more uniform frequency sensitivity at the ex¬ 
pense of independent frequency bins, recovering Fourier phase information which is normally 
lost when constructing the power spectrum from the Fourier amplitudes. 


Figure 4 shows a 5-times overresolved power spectrum of 4U 1626-67 in the immediate 
vicinity of the main pulsation. The power spectrum has been normalized relative to the 
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mean noise power. To improve the signal-to-noise ratio, the power spectrum shown is the 
average of power spectra from three separate 1.8-hr observations (runs 15U, 17B, and 23R). 
In addition to the fundamental pulsation at 130.4398(37) mHz, several other features are 
present. In considering these features, it is important to recall that the peaks due to periodic 
signals in a hnite-length power spectrum are always accompanied by sidelobes. This can be 
understood by considering a hnite time series as the product of an inhnite data stream and a 
rectangular (boxcar) observation window function (equal to unity during the observation and 
zero elsewhere). By the Fourier convolution theorem, the Fourier transform of this product 
is equal to the convolution of the Fourier transforms of the inhnite data stream and the 
observation window function (e.g.. Press et al. 1992). Thus the power at frequency v due to 
a sinusoidal signal with frequency z/q in an ungapped observation of length T is 


P{v) = Po 


simr (uq — iy)T ^ 
7r(z/o - u)T 


( 4 ) 


where Pq is the power at the signal frequency. The dotted curve in Figure 4 shows the ht 
of this function to the data (with the expected mean noise level of unity added), assuming 
that z/q = 130.4398 mHz corresponds to the main pulsation peak. The expected sidelobe 
structure almost exactly accounts for the two symmetric sidelobes on either side of the main 
peak. 


However, after accounting for the expected sidelobe structure of the main pulsation, 
a residual power excess is clearly evident at z/ = 130.045(14) mHz, downshifted 0.395(15) 
mHz from the main pulsation peak. This is consistent with the 0.4011(21) mHz shift for 
the lower sideband reported by Middleditch (1981). The amplitude of the sideband relative 
to the main pulsation is 25% (compared to 20% measured by Middleditch et al. 1981), 
corresponding to a sideband pulsed fraction of 1.5%. To evaluate the statistical signihcance 
of this excess power, note that the powers in an average of three background-dominated 
power spectra obey (within a constant scaling factor) a y^-distribution with six degrees of 
freedom. For a power spectrum normalized such that the mean noise power is unity, the 
probability of a noise fluctuation in a single bin exceeding a threshold power Pthresh is 

Pr(P > Pthresh) = - / (5) 

lo ^ePthresh 

Thus, the observed excess power of 5.9 has a probability of 3.6 x 10“® of being due to a 
random fluctuation, making it significant at the 4.6(T level. Accounting for the fact that all 
bins in the range 0.129-0.132 Hz were searched reduces the statistical significance to d.Ou. 
It is especially striking to note that the frequency shift of the sideband agrees with that 
reported by Middleditch et al. (1981) despite the fact that the actual value of the main 
pulsation frequency (corresponding to the pulsar’s spin) has changed from 130.26 mHz to 



130.44 mHz in the 15 years between the two observations (Chakrabarty et ah 1997). Our 
sideband detection strongly corroborates the 42-min binary period reported by Middleditch 
et al. (1981). Unfortunately, the signal-to-noise ratio of this detection is insufficient to 
permit a more detailed investigation of the binary parameters. 


4. DISCUSSION 

4.1. Multiwavelength timing signatures as probes of LMXBs 

QPOs in accretion-powered pulsars are widely believed to arise from interactions be¬ 
tween the pulsar magnetosphere and the inner edge of an accretion disk. In the beat fre¬ 
quency model (Alpar & Shaham 1985; Lamb et al. 1985), the QPO arises when inhomoge¬ 
neous “clumps” of material at the inner edge of the disk are captured by the magnetosphere 
and accreted by the pulsar, producing a broad power spectral feature at the beat frequency 
^'QPO = — ^^spin between the Keplerian frequency at the inner disk edge and the pulsar’s 

spin frequency. [An alternative explanation, the Keplerian frequency model of van der Klis 
et al. (1987), predicts i'qpo = ^k- However, it can only explain QPOs where zzqpo > z/gpin, 
and thus is not applicable in 4U 1626-67.] Thus, the 0.048 Hz X-ray QPO in 4U 1626-67 
bolsters the premise that this is a disk-fed binary. 

The discovery of an optical QPO at the same frequency shows that a signihcant fraction 
of the optical luminosity from the binary simply mirrors the X-ray emission, and must there¬ 
fore result from X-ray reprocessing. As noted by Middleditch et al. (1981), the existence of 
two distinct, persistent, coherent optical pulsation frequencies suggests that the reprocessing 
occurs in two separate sites, most likely the surface of the mass donor and the accretion 
disk. Simultaneous X-ray and optical timing of these pulsations and the QPO may therefore 
provide a powerful probe of the geometry of the binary and physics of the accretion disk 
(e.g., Arons & King 1993). 

Unfortunately, the BATSE hard X-ray measurements acquired simultaneously with the 
optical observations require 1 d to detect the X-ray pulsations, making them useless for 
a correlative study on short time scales. However, more sensitive X-ray timing missions 
(e.g. RXTE) are well suited to such a task if the observations can be coordinated with 
ground-based optical observations. A comparison of Figures 2 and 3 gives a flavor of the 
kind of data sets that can be directly intercompared if measured simultaneously. Ultravi¬ 
olet timing observations with the Hubble Space Telescope may also provide an important 
probe, as the reprocessing in the disk should be more efficient at these wavelengths than 
in the optical (Arons & King 1993; Anderson et al. 1997). More generally, time-resolved 
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ultraviolet/optical/infrared photometry may provide an important test for candidate multi¬ 
wavelength counterparts in time-variable LMXBs whose companions remain unidentihed. 


4.2. Nature of the mass donor and the mass transfer 

While X-ray binaries and cataclysmic variables with a hydrogen-rich mass donor cannot 
reach orbital periods Porb ^ 80 min (Paczynski & Sienkiewicz 1981), those with a hydrogen- 
depleted mass donor can evolve to extraordinarily short (Porb ~10 min) orbital periods while 
maintaining high mass transfer rates (Nelson, Rappaport, & Joss 1986). Few such systems 
are known, including the six AM CVn cataclysmic variables (see Warner 1995) and the X- 
ray bursters 4U 1820-30 (Porb = 11 min) and 4U 1916-05 (Porb = 50 min). The recently 
reported 20.6-min ultraviolet modulation in the X-ray burster 4U 1850-087 is probably due 
to the same phenomenon (Homer et ah 1996). With a 42-min binary period, 4U 1626-67 is 
clearly a member of this class as well (Levine et al. 1988). 

In addition to the inferred 42-min binary period, there are several other constraints on 
scenarios for the mass donor in 4U 1626-67. The stringent X-ray timing limit (a^jsini < 8 
It-ms; Levine et al. 1988; Shinoda et al. 1990) on the projected radius of the pulsar orbit 
implies 

sin i < 7.8 X 10“^ (1 + , ( 6 ) 

where q = M^/M^ is the mass ratio, M* is the companion mass, M 1.4 is the neutron star mass 
Mx in units of 1.4Mq, and P 42 is the binary period in units of 42 min. The a priori probability 
of observing a system within the inclination limit i < Vax is given by (1 — cos Vax)- 

A second constraint is the assumption that the mass transfer in an ultracompact binary 
is driven by angular momentum losses via gravitational radiation (e.g., Levine et al. 1988). 
The resulting mass transfer rate will be (e.g., Verbunt & van den Heuvel 1995 and references 
therein) 

M = 8.3 X 10-s Ml'! (1 + ^)-i/3 g + | (7) 

where n is the power-law exponent in the companion’s mass-radius relation (P* oc M”), and 
conservative mass transfer via Roche lobe overflow has been assumed. Since the bolometric 
X-ray flux from 4U 1626-67 has been measured (Pj, = 2.4 x 10“® erg cm“^ s“^ in the 0.5-60 
keV band; Pravdo et al. 1979), this also constrains the source distance and X-ray luminosity 
Lx = GMxM/Rx = dTrP^Pj;, where Rx ~ 10 km is the neutron star radius and D is the 
distance to the source. 

Another constraint is a lower limit on M from the observed spin history of the pulsar 
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(Chakrabarty et al. 1997). If the pulsar is spinning up, then the maximum possible accretion 
torque on the neutron star will occur when the magnetospheric radius (where the magnetic 
energy density is comparable to the kinetic energy density of the accreting matter) is close to 
the corotation radius Tco (where the magnetic held lines move with the local disk velocity). 
This limit may be written as N ^ M^/GM^rco, where Vco = (G'Ma,/ 47 r^z/^)^/^ for a Keplerian 
disk. However, the torque on the neutron star can be rewritten as N = 271 where 
Jj. {2/5)MxR‘^ is the pulsar’s moment of inertia (Ravenhall Sz Pethick 1994) and z> is the 
spin frequency derivative - a directly measurable quantity. Noting that z> = 8.5 x 10“^^ Hz 
s“^ was measured during the spin-up of 4U 1626-67, Chakrabarty et al. (1997) thus inferred 
that M ^2x 10"^° Mq yr'h 

Given the measured X-ray hux, the M limit requires D ^3 kpc. In view of the source’s 
high Galactic latitude (6 = —13°), the resulting distance out of the plane is also a considera¬ 
tion when contemplating mass donors. However, this needn’t be a strong consideration, since 
it is plausible for neutron star LMXBs to have high velocities due to a supernova “kick”, 
leading to a considerable distance from the plane over the long X-ray lifetime of the source. 
Indeed, the high-latitude (6 = 38°) LMXB pulsar Her X-l/HZ Her is 6.6 ± 0.4 kpc distant, 
over 4 kpc out of the plane (Reynolds et al. 1997). 

The possible types of Roche-lobe-hlling mass donor for 4U 1626-67 have been discussed 
previously by Levine et al. (1988) and Verbunt, Wijers, & Burm (1990). I review and update 
those discussions here. There are three possible types of Roche-lobe-hlling donors for a 42- 
min neutron star binary. The hrst possibility is a 0.02Mq helium or carbon-oxygen white 
dwarf. From equation (^), the binary inclination for this case must be z ^ 33°, with an a 

_1 /Q 

priori observation probability of 16%. The mass-radius relation for such a star is R* x M* ' 
(Shapiro & Teukolsky 1983), so the mass transfer rate from equation (J^) is 3 x 10 “^^Mq yr“^. 
This corresponds to a distance of about 1 kpc. Although the this case has the most favorable 
inclination constraint of the Roche-lobe-hlling options (due to the extremely small M*), the 
mass transfer rate driven by gravitational radiation alone is an order of magnitude below 
the minimum mass accretion rate inferred by Chakrabarty et al. (1997). 

The second Roche-lobe-hlling possibility is a 0.08Mq hydrogen-depleted star which is 
not fully degenerate, with z ^ 8 ° and an a priori probability of 1%. Assuming R* x M*, 
equation (J^) gives M = 2 x 10“^°Mq yr“^, corresponding to a distance of 3 kpc. These are 
identical to the lower limits on M and D inferred from the pulsar spin-up measurements, and 
thus are just consistent with the observations. Finally, a O. 6 M 0 Roche-lobe-hlling helium¬ 
burning star is also a possible donor for a 42-min neutron star binary. However, the small 
allowed inclination range (z < 1.3°) has an a priori observation probability of only 0.03%, 
and the expected mass transfer rate (M = 3 x IO^^Mq yr“^; Savonije, de Kool, & van den 
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Heuvel 1986) would require an unphysically large distance (36 kpc) to be consistent with the 
observed X-ray flux. 

Angelini et al. (1995) have proposed an alternative to a Roche-lobe-fllling donor. These 
authors detected an unusually strong complex of neon emission lines in the X-ray spectrum 
of 4U 1626-67. Since neon is a product of helium burning, Angelini et al. (1995) suggest 
that the mass donor is a helium-burning star which underfills its Roche lobe, transferring 
matter to the neutron star via a radiatively-driven wind. This possibility is of particular 
interest given the suggestion that the current spin-down of the pulsar is due to accretion 
from a retrograde disk (Nelson et al. 1997), since formation of such a disk is much more 
plausible for wind accretion than for Roche lobe overflow. 

However, the wind from a low-mass helium-burning star would not provide a high enough 
mass transfer rate to be consistent with the X-ray data, as is now shown. A mass of 0.3Mq 
is required for core helium burning (Kippenhahn & Weigert 1990). A O.3M0 helium-burning 
donor would All about 25% of its Roche lobe (Eggleton 1983) and would have R* = 0.077?© 
and an effective temperature T* = 2 x 10^ K (Savonije et al. 1986), yielding a blackbody 
spectrum consistent with or fainter than the optical photometry for distances 4 kpc. More 
generally, helium-burning donors with M* ^ O.bM© would underfill the Roche lobe and be 
consistent with or fainter than the optical measurements for somewhat larger distances, with 
7?* ^ O.147?0 and T* ^ 4x 10^ K. Assuming a radiatively-driven wind similar to that observed 
in other hot stars, the wind mass loss rate with these parameters would be —Mu, ^ lO^^^M© 
yr“^ (Abbott 1982). Especially in light of the poor accretion efficiency expected for a wind- 
fed system (< 0.01; e.g., Frank, King, & Raine 1992), this is several orders of magnitude too 
small to be consistent with the minimum accretion rate, M < 2 x lO^^^M© yr“^. Similar 
considerations seem to rule out any plausible donor transferring matter via a conventional 
radiatively-driven wind. On the other hand, a self-excited wind arising from X-ray heating 
of the mass donor might be a significant contributor to the mass transfer rate, whether or 
not the donor fills its Roche lobe (see Tavani & London 1993 and references therein). 


4.3. Origin of the optical emission 


It is instructive to consider whether a stellar spectrum can explain the optical photom¬ 
etry. The data can in fact be well fit by a single-temperature (stellar) blackbody spectrum, 

27rhz/^/c^ 

772 / exp(/iz//A;T*) — 1’ 



( 8 ) 


where is the flux density from the star per unit frequency z/, 7?* is the stellar radius, D 
is the source distance, h = 6.6 x 10“^'^ erg s is Planck’s constant, c is the velocity of light. 
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k = 1.38 X 10“^® erg K“^ is Boltzmann’s constant, and T* is the effective temperature of the 
star. The best-fit parameters are T* = 1.8 x 10^ K and R^/D = 7.3 x 10“^^, so that D = 
32{R^,/Rq) kpc. A hydrogen-burning star with the correct temperature (a main sequence 
or supergiant B star) is obviously excluded both by the absurd {D ^ 100 kpc) distance 
required and the infinitesimal a priori probability for the necessary binary inclination. A 
O. 3 M 0 He-MS star with D k, 2 kpc would have the correct temperature and flux. However, 
as noted above, it would severely underfill its Roche lobe and so probably could not transfer 
sufficient mass to the pulsar to explain the observed X-ray flux. In addition, the implied 
binary inclination for a 42-min orbit has an a priori probability of < 0.1%. 

However, assuming a typical white dwarf mass-radius relation, {R^/Rq) ^ O.O1(M*/M0)’ 
(Shapiro & Teukolsky 1983), we find D k, 1 kpc for a O.O 2 M 0 white dwarf. As noted above, 
such a star would just fill its Roche lobe, and the expected mass transfer rate would be 
consistent with the observed X-ray flux. For such a close binary. X-ray heating would cer¬ 
tainly dominate the intrinsic luminosity of the degenerate dwarf. Assuming isotropic X-ray 
emission, we have 


'IRI 

1/4 


4 a? TraRl 




= 2.1 X 10" P 42 '^'/lkpc(l + - h*)'/" K, (9) 


where a is the binary separation, rj^ is the X-ray albedo of the white dwarf, and = 
2.4 X 10“® erg cm“" s“" has been used for the X-ray flux (Pravdo et ah 1979). Thus, X-ray 
heating of a low mass WD would naively seem to provide a plausible model for the observed 
optical emission. 


Unfortunately, this is an unrealistic model. The 0.04 Hz QPO detected in the X-ray and 
optical data underscores the presence of an accretion disk (e.g., van der Klis 1995) which 
must contribute significant optical luminosity. The flux from a geometrically thin, optically 
thick disk can be written as 


Fud = 




cos i 


r dr 




exp[hv/kT{r)] — 1 ’ 


( 10 ) 


where T{r) is the disk’s surface temperature as a function of mid-plane radius r. Given the 
X-ray timing limits, it is clear that cosi ~ 1. The outer disk may be assumed to cut off 
sharply at the tidal radius of the neutron star 0.9i?Roche; Frank et al. 1992), corresponding 
Tout ~ 2 X 10 "° cm for 4U 1626-67. The inner accretion disk should be disrupted by the 
pulsar’s magnetosphere at roughly the neutron star’s corotation radius, rjn ^ r^o = 6.5 x 10® 
cm for a Keplerian disk. (Strictly speaking, r^o is an upper limit on rin, but the calculation 
is not sensitive to this quantity since most of the disk’s optical emission comes from its outer 


1/3 
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parts.) For an X-ray heated disk, the disk may be considered as consisting of up to three 
distinct regions: an innermost region powered primarily by internal viscous dissipation, a 
middle region of “shallow” X-ray heating, and an outer region of “deep” X-ray heating (e.g., 
Cunningham 1976; Arons & King 1993). Let us consider each of these regions in turn. 

The standard temperature prohle for an unirradiated thin accretion disk is set by internal 
viscous dissipation (Shakura & Sunyaev 1973; Frank et al. 1992) 

r„ = = 7100 mK* Aft/o rfoi'- K, (11) 

where a = 5.67x 10“^ erg cm“^ K~^ s“^ is the Stefan-Boltzmann constant. Mi 4 is the neutron 
star mass in units of 1.4 Mq, M_iq is the mass accretion rate in units of 10“^° Mq yr“^, and 
rio is the mid-plane radius in units of 10 ^° cm. At sufficiently small radii, this internal heating 
will dominate over X-ray heating in setting the disk temperature. Beyond a critical radius, 
however. X-ray heating of the disk surface will modify the temperature prohle according to 

= + ( 12 ) 


Assuming that the disk is irradiated by a central X-ray point source, then we can write 
as 


T. = 


.hx(l ^d) 
dTTar^ 


COS'lp 


1/4 

'L,{l-Vd)fdH HV 


\dr r j 


-,1/4 


(13) 


where rjd is the X-ray albedo of the disk, ifj is the angle between the normal to the disk 
surface and the vector from the neutron star, and H is the disk’s scale height. Thus, we see 
that the temperature prohle due to X-ray heating is sensitive to the functional form of H (r) 
through the approximation for cosip (which assumes H <^r). 


In the shallow X-ray heating region, the surface temperature prohle is set by X-ray 
heating but the central temperature of the disk is still set by internal viscous dissipation. 
Thus, the disk thickness has the usual value for a standard unirradiated thin disk (Shakura 
& Sunyaev 1973; Frank et al. 1992), 

H = 1.2 X 10® MCf/® m!(o° cm, (14) 


where a is an order unity dimensionless parametrization of the kinematic viscosity u^isc = 
aCgH, Cs is the isothermal sound speed, and jj, is the mean molecular weight in units of the 
hydrogen atomic mass. (Note that /i = 1 for neutral hydrogen, 1/2 for ionized hydrogen, 4 
for neutral helium, 2 for singly-ionized helium, and 4/3 for doubly-ionized helium.) We thus 
hnd 

= 1.2 X lO"* 0-'-!“' (1 - mIP K. 


(15) 
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Shallow X-ray heating will dominate internal viscous heating in setting the disk surface 
temperature when ^ Tq, which occurs at radii exceeding roughly 10® cm. 

If tTq-, then X-ray heating will also dominate internal heating in setting the central 
temperature of the disk (Lyutyi & Sunyaev 1976; Spruit 1995). Here, r is the disk’s optical 
depth, given by r = 34 ^-lo for a standard a-disk (Frank et ah 1992). In this 

region of deep X-ray heating, the vertical structure of the disk will not obey the H oc r®/® 
relationship in equation (pXD, but will instead follow a 77 oc r®/^ power law (e.g., Cunningham 
1976; Vrtilek et ah 1990; Arons & King 1993). However, deep X-ray heating will dominate 
over shallow heating only for radii exceeding roughly 2 x 10^® cm. In a binary as compact 
as 4U 1626-67, the disk will terminate before deep heating becomes dominant, so that a 
combination of internal viscous heating and shallow X-ray heating [equation (12) along with 
equations (11) and (15)] should be valid over the entire disk. 

Accounting for the presence of an X-ray heated accretion disk, it is clear that an X-ray 
heated low-mass white dwarf at a distance of 1 kpc is untenable for 4U 1626-67, since the 
expected optical emission from the disk would be considerably brighter than the measured 
fluxes. Instead, it is more appropriate to consider the opposite case of the disk providing 
most of the optical luminosity by htting equations (10) and (12) along with equations (11) 
and (15) to the measured optical photometry. The free parameters for the £t are D (or, 
equivalently, L^) and rj^. As is evident from Figure 5, the X-ray heated disk model can 
easily fit the data. 

However, as shown in Figure 6, the optical data do not provide a strong joint constraint 
on D and r]d- Ultraviolet observations will eventually provide a much stronger constraint, 
since the UV spectrum from an X-ray heated disk is much more sensitive to the values of 
and r]^- In the absence of ultraviolet observations, it should be noted that de Jong, van 
Paradijs, & Augusteijn (1996) have concluded that LMXB accretion disks have a very high 
effective X-ray albedo (p^ 0.9), based on observations of 4U 1254-69, 4U 1755-33, and 

Sco X-1. If this is the case, then the optical data constrain the distance to 5 ^ ^ 13 kpc 

(95% confidence), corresponding to an X-ray luminosity of 7 x 10®® ^ ^ 5 x 10®^ erg s“^. 
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Table 1. Long Observations of 4U 1626-67 


Run ID 

Date 

Start time 
(UT) 

Duration 

(hr) 

Filter 

High-frequency 

contaminant?’^ 

Pulsations? 

15U 

1995 May 26 

0408 

2.0 

U 

No 

Yes (strong detection) 

17B 

1995 May 26 

0638 

1.8 

B 

Yes 

Yes (strong detection) 

211 

1995 May 27 

0018 

1.7 

I 

No 

No (clouds) 

23R 

1995 May 27 

0308 

1.9 

R 

No 

Yes (strong detection) 

25V2 

1995 May 27 

0604 

0.6 

V 

Yes 

No (clouds) 

29U 

1995 May 28 

0001 

0.7 

U 

Yes 

No (clouds) 

SOU 

1995 May 28 

0046 

1.0 

U 

Yes 

Yes 

32V 

1995 May 28 

0224 

1.0 

V 

Yes 

Yes 

35U 

1995 May 28 

0709 

1.0 

U 

Yes 

Yes 

361 

1995 May 28 

0813 

1.0 

I 

Yes 

Yes (weak detection) 


'^Instrumental signals at harmonics of 1.75 Hz and 60 Hz 


Table 2. Optical Photometry of 4U 1626-67 


Quantity 

U 

B 

V 

R 

/ 

Source count rate (count s“^) 

137(2) 

122 (2) 

111(4) 

177(3) 


Pulsed count rate (count s“^) 

7.9(3) 

6.9(3) 

7.0(5) 

11.8(5) 

25(2) 

Sky count rate (count s“Q 

168(2) 

182(2) 

317(4) 

785(3) 

5405(30) 

Magnitude 

17.50(29) 

18.70(24) 

18.68(15) 

18.68(26) 


Pulsed fraction 

5.8(2)% 

5.7(3)% 

6.3(5)% 

6.7(3)% 


QPO fractional RMS strength 

3.4% 

2.7% 


1.9% 


QPO centroid frequency (Hz) 

0.048 

0.049 


0.049 


QPO FWHM (Hz) 

0.013 

0.012 


0.009 
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FIGURE CAPTIONS 

FIGURE 1: Optical pulse profiles as a function of wavelength for 4U 1626-67. Also shown 
is the 20-60 keV X-ray pulse prohle measured contemporaneously by BATSE. Two pulses 
are shown in each panel, and the pulse phase is plotted relative to the intensity minimum of 
the X-ray pulse. For the /-band optical and BATSE X-ray measurements, the intensity of 
the unpulsed (DC) component of the source emission was not measured. 

FIGURE 2: Fourier power density spectrum of a 2-hour U-band observation (run 15U) of 4U 
1626-67, normalized relative to the mean source power. The dotted line indicates the noise 
power level due to Poisson counting statistics, which has been subtracted off. A 0.048 Hz 
quasi-periodic oscillation with 4.9% fractional RMS amplitude is clearly detected, as is the 
7.66-s coherent pulsation due to the neutron star’s rotation. An underlying 1/z/ spectrum is 
measured on time scales ^ 10 s. 

FIGURE 3: Fourier power density spectrum of a 71 ks 2-6 keV RKTE/VCA observation 
of 4U 1626-67 on 1996 Feb 11, normalized relative to the mean X-ray source power. The 
dotted line indicates the noise power level due to Poisson counting statistics, which has 
been subtracted off. Although this power spectrum cannot be directly compared to the 
optical power spectrum in Figure 2 (since they were not acquired contemporaneously), the 
overall similarities are obvious and bolster the idea that the optical emission is dominated 
by rapid X-ray reprocessing. Several higher harmonics of the 0.1304 Hz coherent pulsation 
are observed due to the highly non-sinusoidal X-ray pulse shape in the 2-6 keV band. 

FIGURE 4: A 5x over-resolved Fourier power spectrum of 4U 1626-67 in the vicinity of 
the coherent 7.66-s pulsation, normalized relative to the local noise power. The solid line 
is the average of three separate power spectra of 1.8-hr observations. The dotted line is 
a model of the main pulsation’s expected sidelobe due to the hnite length of the power 
spectrum. This model is able to explain the symmetric sidelobes on either side of the main 
pulsation. A residual sideband downshifted 0.395(15) mHz from the main pulsation (25% 
relative amplitude) is also present, with a statistical signihcance of 4(T. This frequency shift 
is consistent with a 2530 ± 100 s prograde binary orbit. 

FIGURE 5: Optical {UBVRT) photometry of 4U 1626-67 on 1995 May 26. The solid curve 
indicates a typical X-ray heated disk solution. The particular solution shown assumed 7]^ = 
0.95, D = 8.5 kpc, = 6.5 x 10® cm, and rout = 2 x 10^° cm. However, a range of values 
for rjd and D provide a good £t to the data (see Figure 6). 

FIGURE 6: Distance to 4U 1626-67 as a function of effective X-ray albedo of the disk, 
assuming that all the optical flux is due to an X-ray heated accretion disk with rjn = 6.5 x 10® 
cm, and rout = 2 x 10^'^ cm. The shaded region is consistent with the optical photometric 
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data at the 95%-confidence level. Based on observations of other LMXBs, de Jong et ah 
(1996) conclude that rjd ^ 0.9. If this is correct, then the solutions allowed by the data must 
lie in the shaded region to the right of the vertical dotted line. 
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